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A  series  of  Li  deficient  LiTi2(P04)3  samples  were  prepared  and  sintered  and  the  density  was  measured  to 
determine  the  rate-controlling  species  for  sintering  of  LiTi2(P04)3.  It  was  observed  that  as  the  Li  content 
decreased  the  density  decreased.  This  result  suggests  that  oxygen  does  not  control  sintering.  A  compar¬ 
ison  of  the  LiTi2(P04)3  sintering  data  to  sintering  and  diffusion  data  in  olivine,  which  exhibits  a  similar 
framework  structure  to  LiTi2(P04)3,  suggests  that  P  is  the  species  which  controls  sintering.  This  suggestion 
was  confirmed  by  the  density  results  of  a  Li  excess  LiTi2(P04)3  sample. 
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1.  Introduction 

Recently,  there  has  been  a  renewed  interest  in  the  develop¬ 
ment  of  high  energy  Li-air  batteries.  One  configuration  involves 
the  use  of  a  Li  anode  in  a  non-aqueous  electrolyte  which  is  sep¬ 
arated  from  an  aqueous  electrolyte  containing  the  air  cathode  by 
a  solid-state  Li-ion  conducting  membrane  [1].  One  of  the  require¬ 
ments  for  the  membrane  material  besides  high  Li-ion  conductivity 
is  that  it  must  also  be  chemically  stable  with  water;  as  a  result  Li- 
ion  conductors  based  on  an  oxygen  framework  are  the  preferred 
choice.  In  this  case,  there  are  essentially  three  possibilities,  those 
based  on  a  pervoskite  [2,3],  garnet  [4,5]  or  Na  super  ion  conduc¬ 
tor  (NASICON)  [6-11  ]  structure.  Of  these  choices,  LiTi2(P04)3  (LTP), 
based  on  the  NASICON  structure,  is  one  of  the  most  widely  inves¬ 
tigated  solid-state  Li-ion  conducting  materials  [6-11].  For  LTP  it 
has  been  observed  that  Li-ion  conductivity  is  a  strong  function 
of  density.  For  example,  Li-ion  conductivity  can  vary  from  about 
2  x  10_6Scm_1  [9]  for  a  material  with  a  relative  density  (actual 
density/theoretical  density)  around  70%  to  1  x  10-3  S  cm-1  for  a 
material  with  a  relative  density  close  to  100%  [13].  It  has  been 
suggested  that  high  density  is  required  to  reduce  grain  boundary 
resistance  [6].  In  order  to  achieve  the  high  densities  two  process¬ 
ing  approaches  are  commonly  used.  The  first  approach  is  melting 
of  doped-LTP  followed  by  quenching  to  form  a  glass  which  is  sub¬ 
sequently  reheated  to  allow  for  the  formation  of  a  glass-ceramic 
[  12,13  ].  The  second  approach  is  to  use  conventional  sintering  [6-11  ]. 
It  has  been  observed  that  pure  LiTi2(P04)3  cannot  be  sintered  to 
high  densities  [6-11,14,15].  In  order  to  obtain  a  dense  material 
LiTi2(P04)3  must  be  doped  with  a  trivalent  impurity  on  the  Ti  site  or 
a  low  melting  phase  added  [6-11  ].  In  both  of  these  cases  high  den¬ 
sity  is  achieved  as  a  result  of  a  second  (i.e.,  A1P04)  or  liquid  phase 
along  the  grain  boundaries,  which  enhances  sintering  [6-11,14,15]. 
However,  at  present  there  is  no  data  or  explanation  in  the  litera¬ 
ture  on  what  controls  the  sintering  behavior  of  pure  LiTi2(P04)3 
and  hence,  why  a  high  density  material  cannot  be  obtained  after 
sintering. 


For  the  case  of  sintering  it  is  controlled  by  the  concentration  and 
mobility  of  the  slowest  species  [16-19].  For  the  case  of  most  oxides 
it  is  the  oxygen  ion  which  is  the  slowest  moving  species  and  moves 
via  a  vacancy  mechanism  and  hence,  controls  sintering  [16-20]. 
Thus,  for  LiTi2(P04)3  it  is  expected  that  oxygen  vacancies  will  con¬ 
trol  the  sintering  behavior.  One  possible  way  to  test  this  hypothesis 
is  to  alter  the  concentration  of  oxygen  vacancies  in  LiTi2(P04)3  and 
see  how  this  effects  sintering.  For  example,  when  Li  vacancies  [V'J 
(Kroger-Vink  [16-19]  notation  is  used)  are  created  in  LiTi2(P04)3, 
then  according  to  the  charge  neutrality  condition  for  oxidizing  con¬ 
ditions  (i.e.,  heat-treatment  under  an  air  atmosphere)  given  in  Eq. 
(1 )  an  increase  in  oxygen  vacancies  [  V0**  ],  will  result  and  hence,  an 
increase  in  the  sintering  rate  and  a  higher  density  will  be  obtained. 

Ki]  =  2[\70"]  (1) 

It  is  the  purpose  of  this  note  is  to  determine  the  rate-controlling 
species  for  sintering  of  LiTi2(P04)3  by  varying  the  Li  vacancy  con¬ 
centration  (i.e.,  oxygen  vacancy  concentration)  and  observing  its 
affect  on  density. 

2.  Experimental 

Different  Li  vacancy  concentration  LTP  samples  were  prepared 
by  using  Li  deficient  LTP.  Li!_xTi2(P03)4  samples  with  x  =  0.0,  0.05, 
0.10,  and  0.15  samples  were  chosen.  The  exact  formulae  for  the 
lithium  deficient  LTP  materials  is  of  the  form;  Lii_xTi2P3012_x/2, 
which  incorporates  the  charge  compensating  oxygen  vacancies. 
However,  for  the  rest  of  the  paper  the  different  sample  composi¬ 
tions  will  be  represented  by  the  following  formulas,  LiTi2(P04)3, 
Lio.95Ti2(P04)3,  Li090Ti2(PO4)3  and  Li0.85Ti2(PO4)3.  This  does  not 
affect  the  results  and  conclusions.  All  samples  were  synthesized 
using  a  solid-state  method  from  Ti02  (rutile  structure),  Li2C03  and 
NH4H2P04.  The  starting  materials  were  ball  milled  using  methanol 
for  6  h.  The  powders  were  ground  into  a  fine  powder  and  calcined 
at  900  °C  under  air  for  4  h  in  a  covered  alumina  crucible.  The  result¬ 
ing  powder  was  reground  and  uniaxially  pressed  into  cylindrical 
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specimens  13  mm  in  diameter  with  a  thickness  between  2  and 
4  mm  at  pressure  of  200  MPa.  A  5  wt.%  polyvinyl  alcohol  binder 
was  used.  The  pellets  were  then  sintered  at  950  °C  for  2  h  under  air 
on  platinum  foil.  The  Li  content  of  the  samples  after  sintering  was 
determined  using  inductively  coupled  plasma  analysis  (Galbraith 
Laboratories,  Inc.,  Knoxville,  TN).  All  samples  were  characterized 
by  X-ray  diffraction.  Lattice  constants  were  determined  by  fitting 
the  data  using  Rietveld  refinement  [21  ].  The  density  of  the  sintered 
pellets  was  determined  from  the  weight  and  physical  dimensions. 

3.  Results  and  discussion 

The  LiTi2(P04)3,  Lio.g5Ti2(P04)3t  Lio.goTi2(P04)3  and 
Lio.85Ti2(P04)3  samples  were  all  white  colored  after  sintering. 
The  X-ray  diffraction  pattern  of  the  LiTi2(P04)3  sample  after 
sintering  is  shown  in  Fig.  1.  From  Fig.  1  it  can  seen  that  there  are 
two  phases  present,  a  majority  LiTi2(P04)3  phase  and  a  TiP207 
second  phase.  Second  phase  TiP207  was  found  in  all  four  samples. 
The  amount  of  TiP207  determined  using  Rietveld  refinement  was 
between  4  and  6  wt.%  for  all  four  samples.  TiP207  is  expected 
according  to  Eq.  (2)  as  a  result  of  Li20  loss  during  heating  and 
has  been  observed  in  many  LTP  samples  after  calcining  and/or 
sintering  [15,22,23].  The  X-ray  diffraction  patterns  did  not  reveal 
the  presence  of  Ti02. 

2LiTi2(P04)3  -*  3TiP207  +Ti02  +  t  Li20  (2) 

Inductively  coupled  plasma  analysis  did  confirm  Li  loss  dur¬ 
ing  heat-treatment.  The  actual  Li  content  of  the  LiTi2(P04)3, 
Lio.95Ti2(P04)3,  Li0.goTi2(P04)3  and  Li0.85Ti2(PO4)3  samples  after 
heat-treatment  was  0.96  ±0.02,  0.90  ±0.02,  0.85  ±0.02  and 
0.81  ±0.02  instead  of  the  initial  values  of  1.0,  0.95,  0.90  and  0.85, 
respectively.  The  fact  that  all  the  samples  had  the  same  Li  loss 
(0.04-0.05)  after  heat-treatment  is  expected  since;  they  were  all 
given  the  exact  same  heat-treatment. 

The  X-ray  diffraction  pattern  for  the  majority  LiTi2(P04)3  phase 
can  be  indexed  to  a  rhombohedral  (R3c)  NASICON  structure  for 
all  four  samples.  The  lattice  parameters  a  and  c,  determined  from 
Rietveld  refinement  as  a  function  of  Li  content  (amount  determined 
from  the  inductively  coupled  plasma  analysis)  are  shown  in  Fig.  2 A 
and  B,  respectively.  From  Fig.  2  several  important  points  are  noted. 
First,  the  lattice  parameters  for  LiTi2(P04)3  are  in  good  agreement 
with  those  in  the  literature  [7,10,23-26].  Secondly,  both  the  a  and 
c  lattice  parameters  decrease  with  decreasing  Li  content  (xt).  This 
result  is  consistent  with  the  formation  of  vacancies  in  the  Li  defi- 


Fig.  2.  Lattice  parameters  of  Lii_xTi2(P04)3  as  a  function  of  x  (using  the  Li  content 
after  sintering). 

cient  LTP  materials.  We  would  expect  that  from  Eq.  (1)  that  with 
decreasing  Li  content  (xt)  an  increase  in  the  number  of  oxygen 
vacancies  and  hence,  an  increase  in  density. 

Fig.  3  is  a  plot  of  relative  density  versus  Li  content  (amount 
determined  from  the  inductively  coupled  plasma  analysis)  for 
the  LiTi2(P04)3,  Lio.gsTi2(P04)3,  Lio.goTi2(P04)3  and  Lio.8sTi2(P04)3 
samples.  The  relative  density  is  the  actual  density  obtained  from  the 
weight  and  physical  dimensions  divided  by  the  theoretical  density 
of  the  sample,  ptheo  ample,  given  in  the  following  equation: 

Ptheo  sample  —  AheoLTp/LTP  +  PtheoTiP207 /liP207  (3) 

where  Aheo  ltp  is  the  theoretical  density  of  LTP  (changing  from 
LiTi2(P04)3,  to  Li0.85Ti2(PO4)3  has  <1%  difference  in  density,  hence 
the  density  value  for  LTP  from  the  JPDS  X-ray  PDF  file,  #00-35- 
754  was  used,  2.948  g  cm-3 ),  ptheoTiP2o7  is  the  theoretical  density  of 
TiP207  (3.014 gem-3,  JPDS  X-ray  PDF  file  #00-038-1468), /LTP  is  the 
fraction  of  LTP  and/Tip2o7  is  the  fraction  of  TiP207.  For  the  four  dif¬ 
ferent  LTP  samples  since,  the  amount  of  TiP207  was  nearly  the  same 
in  all  four  samples  (4-6  wt.%),/np2o7  =  5  wt.%  and /LTP  =  95  wt.%  was 
used  in  calculating  ptheo  sample-  Inserting  these  values  into  Eq.  (3) 
results  in  95%  of  the  sample  density  coming  from  the  LTP  phase 
thus,  any  trend  in  the  density  data  with  Li  content  must  be  due  to 
the  behavior  of  the  LTP  phase.  From  Fig.  3  several  important  points 
are  noted.  First,  the  relative  density  of  LiTi2(P04)3  is  69%.  This  value 
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Fig.  3.  Relative  density  of  Lii_xTi2(P04)3  as  a  function  of  x  (using  the  Li  content  after 
sintering). 

is  in  good  agreement  with  the  value  of  66%  for  LiTi2(P04)3  sin¬ 
tered  at  almost  the  same  temperature  (900  °C)  for  the  same  time 
(2h)  recorded  by  Aono  et  al.  [8].  It  is  also  in  agreement  with  the 
value  of  69%  observed  by  Yamamoto  et  al.  [27]  for  LiGe2(P04)3  sin¬ 
tered  at  900  °C  for  2h.  Secondly,  as  the  Li  content  (xf)  increases 
the  relative  density  decreases.  This  result  reveals  as  the  number 
of  Li  vacancies  increases  and  hence,  the  number  of  oxygen  vacan¬ 
cies  increases  a  decrease  in  density  is  exhibited.  This  is  opposite  of 
what  is  expected  if  oxygen  vacancies  control  sintering.  For  the  case 
of  sintering  controlled  by  oxygen  vacancies  an  increase  in  density 
should  be  observed  as  the  number  of  oxygen  vacancies  increases. 
This  result  suggests  oxygen  vacancies  are  not  the  rate-controlling 
species.  This  observation  is  not  unreasonable  when  the  LTP  sinter¬ 
ing  results  are  compared  to  sintering  and  diffusion  data  for  olivine 
materials.  The  basic  LTP  structure  consists  of  a  three-dimensional 
framework  of  Ti06  octahedra  interconnected  with  P04  tetrahedra 
[6,25,26].  The  basic  olivine  structure  can  be  characterized  by  M06 
octahedra  (where  M  is  a  metal)  interconnected  with  Si04  tetrahe¬ 
dra  [16-18].  Since,  both  LTP  and  olivine  have  a  similar  structure, 
the  sintering  and  diffusion  data  for  olivine  can  be  used  as  a  guide  to 
determine  the  rate-controlling  species  for  the  sintering  of  LTP.  For 
the  case  of  olivine,  diffusion  and  sintering  studies  have  revealed 
that  Si  not  oxygen  is  the  slower  moving  species  and  hence,  con¬ 
trols  sintering  [28-30].  Based  upon  this,  it  is  highly  likely  that  P  is 
the  species  that  controls  the  sintering  rate  of  LTP.  If  this  is  true  and 
assuming  that  the  ionic  defects  in  LTP  form  Schottky  defects  (vacan¬ 
cies)  rather  than  Frenkel  defects  (interstitials),  then  the  Schottky 
defect  equilibria  reaction  for  LiTi2(P04)3  is  given  in  the  following 
equation: 

Kj]  [V"']2[V"']3[Vo-]12  =  ks  (4) 

where  [\7Ti////]  is  a  titanium  vacancy,  [VP/////]  is  a  phosphorous  vacancy 
and  J<s  is  the  Schottky  equilibrium  constant.  From  Eq.  (4)  it  is 
observed  that  the  product  of  the  Li,  Ti,  P  and  O  vacancy  concentra¬ 
tions  is  fixed  at  a  given  temperature.  Thus,  if  one  concentration  is 
changed  the  others  must  vary  to  keep  this  product  constant.  For  our 
case  we  have  increased  both  the  Li  and  O  vacancy  concentrations  as 
LTP  becomes  more  Li  deficient,  as  a  consequence  according  to  Eq. 
(4)  the  Ti  and  P  vacancy  concentrations  must  decrease.  As  a  result  of 
the  decreased  P  vacancy  concentration,  assuming  that  P  is  rate  con¬ 
trolling,  we  will  have  lower  sintering  rate  and  hence,  a  lower  density 
will  be  achieved.  The  results  shown  in  Fig.  3  reveal  a  lower  density  is 


achieved  as  the  Li  and  O  vacancy  concentrations  increase,  in  agree¬ 
ment  with  the  suggestion  that  P  is  the  rate-controlling  species  for 
the  sintering  of  LTP. 

To  confirm  this  suggestion,  a  LTP  sample  with  excess  Li  was 
prepared.  A  Li  excess  Li1+xTi2(P03)4  sample  with  a  starting  x  =  0.10 
was  prepared  and  sintered  under  the  same  conditions  as  that  for 
the  LiTi2(P04)3,  Li0.95Ti2(PO4)3,  Li0.goTi2(P04)3  and  Li0.8sTi2(PO4)3 
samples.  X-ray  diffraction  of  the  Li! 4Ti2(P04)3  sample  also  revealed 
not  only  the  majority  LiTi2(P04)3  phase  but  also  a  TiP2Oy  second 
phase,  just  like  for  the  case  of  the  Li  deficient  LTP  samples.  Rietveld 
refinement  revealed  the  amount  of  TiP207  second  phase  was  4  wt.% 
in  the  Lii.iTi2(P04)3  sample,  an  amount  similar  to  that  observed  for 
the  Li  deficient  LTP  samples.  The  final  Li  content  of  the  Li! 4Ti2(P04)3 
sample  was  not  analyzed.  However,  based  on  the  inductively  cou¬ 
pled  plasma  results  for  Li  deficient  LTP  samples  which  had  a  Li  loss 
(0.04-0.05)  after  heat-treatment  it  is  expected  that  actual  Li  con¬ 
tent  of  the  Li!.!Ti2(P04)3  sample  is  about  1.05  instead  of  1.10,  which 
still  makes  this  sample  Li  excess. 

For  the  case  of  Li  excess  LTP  the  charge  neutrality  condition  is 
given  by  the  following  equation: 

[V]  =  WH  +  5[Vpw/]  (5) 

where  [Li#]  is  the  concentration  of  extra  Li  on  interstitial  sites.  At 
present  it  is  not  possible  to  determine  if  one  type  of  vacancy,  Ti 
or  P  is  dominant.  However,  from  Eq.  (5)  it  can  be  observed  that 
the  extra  Li  results  in  a  increase  in  the  number  of  Ti  and  P  vacan¬ 
cies  in  the  Li  excess  LTP  compared  to  LiTi2(P04)3.  Assuming  P  is 
the  rate-controlling  species  for  sintering  then,  a  higher  density 
is  expected  for  the  Li  excess  LTP  compared  to  LiTi2(P04)3.  The 
Li!  !Ti2(P04)3  sample  had  a  relative  density  of  74%  compared  69% 
for  LiTi2(P04)3.  This  result  adds  further  confirmation  that  P  is  the 
rate-controlling  species  for  sintering  of  LiTi2(P04)3.  The  results  of 
this  study  suggest  that  in  order  to  obtain  highly  dense  LiTi2(P04)3 
without  the  addition  of  secondary  phases  along  grain  bound¬ 
aries  through  conventional  sintering  point  defects  which  increase 
the  concentration  of  P  vacancies,  such  as  adding  excess  Li  are 
needed. 

4.  Conclusions 

A  series  of  Li  deficient  LTP  samples,  LiTi2(P04)3,  Li0.95Ti2(PO4)3i 
Li0.9oTi2(P04)3  and  Li0.85Ti2(PO4)3  were  prepared  and  sintered  and 
the  density  was  measured  to  determine  the  rate-controlling  species 
for  sintering  of  LiTi2(P04)3.  It  was  observed  that  as  the  Li  content 
(i.e.,  the  number  of  Li  vacancies  and  oxygen  vacancies)  decreased 
the  density  decreased.  This  trend  is  opposite  of  what  is  expected 
if  oxygen  vacancies  control  sintering.  A  comparison  of  the  LTP 
sintering  results  to  diffusion  and  sintering  data  in  olivine,  which 
exhibits  a  similar  structure  to  LTP,  suggests  that  P  is  the  species 
which  controls  sintering  of  LiTi2(P04)3.  This  suggestion  was  con¬ 
firmed  by  the  higher  density  of  a  Li  excess  LTP  sample  compared 
to  LiTi2(P04)3.  In  order  to  obtain  highly  dense  LiTi2(P04)3,  without 
the  addition  of  secondary  phases  along  grain  boundaries,  using  con¬ 
ventional  sintering,  point  defects  which  increase  the  concentration 
of  P  vacancies  (i.e.,  for  example,  adding  excess  Li)  are  needed. 

Acknowledgement 

The  authors  would  like  to  acknowledge  support  of  the  U.S.  Army 
Research  Laboratory. 

References 

[1]  S.J.  Visco,  B.D.  Katz,  Y.S.  Nimon,  L.C.  De  Jonghe,  US  Patent  7,282,295  B2  (2007). 

[2]  V.  Thangadurai,  W.  Weppner,  Solid  State  Ionics  12  (2006)  81. 


Letter  to  the  Editor  /  Journal  of  Power  Sources  182  (2008)  626-629 


629 


[3]  Y.  Inaguma,  L.Q,  Chen,  M.  Itoh,  T.  Nakamura,  T.  Uchida,  H.  Ikuta,  M.  Wakihara, 
Solid  State  Commun.  86  (1993)  689. 

[4]  V.  Thangadurai,  S.  Schwenzel,  W.  Weppner,  Ionics  11  (2005)  11. 

[5]  V.  Thangadurai,  W.  Weppner,  J.  Am.  Ceram.  Soc.  88  (2005)  411. 

[6]  G.  Adachi,  N.  Imanaka,  H.  Aono,  Adv.  Mater.  8  (1996)  127. 

[7]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  Solid  State  Ionics  40-41 
(1990)38. 

[8]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  Solid  State  Ionics  47 
(1991)257. 

[9]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  J.  Electrochem.  Soc. 
140(1993)1827. 

[10]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  J.  Electrochem.  Soc. 
137(1990)1023. 

[11]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  Chem.  Lett.  (1990) 
1825. 

[12]  J.  Fu,  US  Patent  6,485,622  B1  (2002). 

[13]  J.  Fu,  Solid  State  Ionics  96  (1997)  195. 

[14]  Y.  Nan,  W.E.  Lee,  P.F.  James,  J.  Am.  Ceram.  Soc.  75  (1992)  1641. 

[  15  ]  S.  Wong,  P.J.  Newman,  A.S.  Best,  K.M.  Nairn,  D.R.  MacFarlane,  M.  Forsyth, J.  Mater. 
Chem.  8  (1998)  2199. 

[  16]  M.V.  Barsum,  Fundamentals  of  Ceramics,  The  McGraw-Hill  Companies  Inc.,  New 
York,  1997. 

[17]  Y.M.  Chiang,  D.  Birnie  III,  W.D.  Kingery,  Physical  Ceramics,  Wiley,  New  York, 
1997. 

[18]  W.D.  Kingery,  H.K.  Bowen,  D.R.  Uhlmann,  Introduction  to  Ceramics,  second  ed., 
Wiley,  New  York,  1976. 

[19]  P.  Kofstad,  Nonstoichiometry,  Diffusion  and  Electrical  Conductivity  in  Binary 
Metal  Oxides,  Wiley-Interscience,  New  York,  1972. 

[20]  H.J.  Frost,  M.F.  Ashby,  Deformation-Mechanism  Maps,  Pergamon  Press,  Oxford, 
1982. 

[21  ]  H.M.  Rietveld,  J.  Appl.  Crystallogr.  2  (1969)  65. 

[22]  Y.  Kobayashi,  T.  Takeuchi,  M.  Tabuchi,  A.  Ado,  H.  Kageyama,  J.  Power  Sources 
81-82(1999)853. 


[23]  N.V.  Kosova,  E.T.  Devyatkina,  A.P.  Stepanov,  A.L.  Buzlukov,  Ionics,  in  press. 

[24]  H.  Aono,  E.  Sugimoto,  Y.  Sadaoka,  N.  Imanaka,  G.  Adachi,  J.  Electrochem.  Soc. 
136(1989)590. 

[25]  V.  Thangadurai,  A.K.  Shukla,  J.  Gopalakrishnan,  J.  Mater.  Chem.  9  (1990) 
739. 

[26]  K.  Arbi,  S.  Mandal,  J.M.  Rojo,  J.  Sanz,  Chem.  Mater.  14  (2002)  1091. 

[27]  H.  Yamamoto,  M.  Tabuchi,  T.  Takeuchi,  H.  Kageyama,  O.  Nakamura,  J.  Power 
Sources  68  (1997)  397. 

[28]  O.  Jaoul,  J.  Geophys.  Res.  95  (17)  (1990)  631. 

[29]  O.  Jaoul,  B.  Houlier,  F.  Abel,  J.  Geophys.  Res.  88  (1983)  613. 

[30]  B.  Houlier,  O.  Jaoul,  F.  Abel,  R.C.  Liebermann,  Phys.  Earth  Planet.  Inter.  50  (1988) 
240. 

J.  Wolfenstine* 
D.  Foster 
J.  Read 
J.L.  Allen 

Army  Research  Laboratory,  AMSRD-ARL-SE-DC, 
2800  Powder  Mill  Road,  Adelphi,  MD  20783,  United  States 

*  Corresponding  author.  Tel:  +1  301  394  0317; 

fax:  +1  301  394  0273. 
E-mail  address:  jwolfenstine@arl.army.mil 
(J.  Wolfenstine) 

31  March  2008 
Available  online  29  April  2008 


